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Abstract 

Extending the notion of indispensable binomials of a toric ideal (0))) we 

define indispensable monomials of a toric ideal and establish some of their properties. 
They are useful for searching indispensable binomials of a toric ideal and for proving 
the existence or non-existence of a unique minimal system of binomials generators 
of a toric ideal. Some examples of indispensable monomials from statistical models 
for contingency tables are given. 

1 Introduction 

In recent years techniques of computational commutative algebra found applications in 
many fields, such as optimization ((H), computational biology (ficl 0; 0), and statistics 



(111). Particularly, the algebraic view of discrete statistical models has been applied in 
many statistical problems, including conditional inference (Q), disclosure limitation (jlih . 
the maximum likelihood estimation (Q), and parametric inference (j^. Algebraic statistics 
is a new field, less than a decade old, and its term was coined by Pistone, Riccomagno 
and Wynn, by the title of their book (fill). Computational algebraic statistics has been 
very actively developed by both algebraists and statisticians since the pioneering work of 
Diaconis and Sturmfels ((3|). For sampling from a finite sample space using Markov chain 
Monte Carlo methods, Diaconis and Sturmfels (j^) defined the notion of Markov bases and 
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showed that a Markov basis is equivalent to a system of binomial generators of a toric 
ideal. 

In statistical applications, the number of indeterminates is often large, and at the same 
time, there exists some inherent symmetry in the toric ideal. For this reason, we encounter 
computational difficulty in applying Grobner bases to statistical problems. In particular, 
even the reduced Grobner basis of a toric ideal may contain more than several thousands 
elements, but one might be able to describe the basis concisely using symmetry. For 
example, i shows that the unique minimal Markov bases for 3x3xi^", > 5 contingency 
tables with fixed two-dimensional marginals contain only 6 orbits with respect to the 
group actions of permuting levels for each axis of contingency tables, while there are 3240, 
12085, and 34790 elements in reduced Grobner bases for K = 5,6, and 7, respectively. 
Furthermore in this example the reduced Grobner basis contains dispensable binomials 
and is not minimal. 

Because of the difficulty mentioned above, the first two authors of this paper have 
been investigating the question of minimality of Markov bases. In (fl3 ). we defined the 
notion of indispensable moves, which belong to every Markov basis. We showed that there 
exists a unique minimal Markov basis if and only if the set of indispensable moves forms 
a Markov basis. Shortly after, Ohsugi and Hibi investigated indispensable binomials. 
They showed that the set of indispensable binomials coincides with the intersection of 
all reduced Grobner basis with respect to lexicographic term orders in (0). Thus, we 
are interested in enumerating indispensable binomials of a given toric ideal. However, in 
general, the enumeration itself is a difficult problem. 

This paper proposes the notion of indispensable monomials and investigate some of 
their properties. The set of indispensable monomials contains all terms of indispensable 
binomials. Therefore if we could enumerate indispensable monomials, then it would be 
straightforward to enumerate indispensable binomials. Here it may seem that we are 
replacing a hard problem by a harder one. Computationally this may well be the case, 
but we believe that the notion of indispensable monomials may be useful for understanding 
indispensable binomials and finding the existence of the unique minimal Markov basis. 

In Section 121 we will set notation and summarize relevant results from (jl4l). In Section 
121 we will define indispensable monomials and prove some basic properties of the indis- 
pensable monomials. Further characterizations of indispensable monomials are given in 
Section HI and some nontrivial examples are given in Section We will conclude with 
some discussions in Section IHl 



2 Preliminaries 



In this section we will set appropriate notation and then summarize main results from (|l4l). 
Because of the fundamental equivalence mentioned in (j^, we use "system of binomial 
generators" and "Markov basis" synonymously. Other pairs of synonyms used in this 
paper are ( "binomial" , "move" ) , ("monomial", "frequency vector" ) and ("indeterminate", 
"cell"), as explained below. 
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2.1 Notation 



Because this paper is based on ()141 ). we follow its notation and terminology in statistical 
context. Also we adapt some notation in (jl^ Ifih. Vectors, through this paper, are column 
vectors and x' denotes the transpose of the vector x. 

Let X be a finite set of p = |X| elements. Each element of X is called a cell. By ordering 
cells, we write X = {1, . . . ,p} hereafter. A nonnegative integer e N = {0, 1, . . .} is a 
frequency of a cell i and x = (xi, . . . , Xp)' G is a frequency vector (nonnegative integer 
vector). We write |x| = Yl^=i denote the sample size of x. In a framework of similar 
tests in statistical theory (see Chapter 4 of we consider a dimensional sufficient 
statistic defined by 

p 

t = ^ ^ ^iXi, 
i=l 

where G Z'^ = {0, ±1, . . Y is a (i- dimensional fixed integral column vector for i = 
1, . . . ,p. Let A = (ai, . . . , ap) = {aji) denote a dx p integral matrix, where aji is the j-th 
element of a^. Then the sufficient statistic t is written as t = ^x. The set of frequency 
vectors for a given sufficient statistic t is called a t- fiber defined by 

= {x G N*' I t = ^x}. 

A frequency vector x (g W) belongs to the fiber JF^^ by definition. We assume that a 
toric ideal is homogeneous, i.e. there exists w such that w'aj = 1, i = 1, . . . ,p. In this 
case the sample size of t is well defined by |t| = |x| where x G J?-t. If the size of is 1, 
i.e. 

^^x = {x}, 

we call X G a 1-element fiber. denotes the size (the number of the elements) 

of the fiber JFax- The support of x is denoted by supp(x) = {i | Xj > 0} and the i-th. 
coordinate vector is denoted by ej = (0, . . . , 0, 1, 0, . . . , 0)', where 1 is in the z-th position. 

Now, we consider the connection between contingency tables and toric ideals. Let 
. . . , Mp] = fc[u] denote the polynomial ring in p indeterminates ui, . . . ,Up over the 
field k. We identify the indeterminate G u with the cell i E I. For a p-dimensional 
column vector x G of non- negative integers, let u'' = Ui^---Up^ G k[u\ denote a 
monomial. For the sufficient statistic t, we also treat t = [ti, . . . ,trf)' as indeterminates. 
Let A;[t^^] = k[ti, . . . ,tci,t^^ , . . . ,t'^^] denote the Laurent polynomial ring. Then the 
system of equations t = ^x is identified as the mapping tt : k[u] fc[t^^] defined as 
Xi H- >■ t*^* = t1^^ ■ ■ -f^/^. The kernel of tt is denoted by = ker(7r) and it is the toric ideal 
associate to the matrix A. 

For statistical applications, it is important to construct a connected Markov chain 
over the given t-fiber. In (j^, Diaconis and Sturmfels showed that a generator of the toric 
ideal Ij^ forms a Markov basis, i.e., it can give a connected chain for any t-fiber. 

A p-dimensional integral column vector z G is a move (for A) if it is in the kernel 
of A, 

Az = 0. 
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Let z+ = {zf , . . . , Zp)' and z~ = {z^ , . . . , z~y denote the positive and negative part of a 
move z given by 

z^ = max{zi, 0), z~ = — min(zj, 0), 

respectively. Then z = z+ — z~ and z+ and z~ are frequency vectors in the same fiber 
^Az+{= ^Az~)- Adding a move z to any frequency vector x does not change its sufficient 
statistic, 

^(x + z) = ^x, 

though x + z may not necessarily be a frequency vector. If adding z to x does not produce 
negative elements, we see that x e JF^^ is moved to another element x + z G JF^^ by z. 
In this case, we say that a move z is applicable to x. z is applicable to x if and only if 
X + z G J-'ax, and equivalently, x > z~, i.e., x — z^ G W. In particular, z is applicable 
to z~. We say that a move z contains a frequency vector x if z+ = x or z^ = x. The 
sample size of z"*" (or z~) is called a degree of z and denoted by 



deg(z) = |z 



+ 1 



We also write |z| = J2i=i k«l — 2deg(z). 

Let B = {zi, . . . , Zi} be a finite set of moves. Let x and y be frequency vectors in the 



same fiber, i.e., x, y G J-'Axi= ^Ay)- Following (|l4f). we say that y is accessible from x by 
B if there exists a sequence of moves Zj^, . . . , Zj^. from B and e^- G {—1, +1}, j = 1, • • • , fc, 
satisfying y = x + Y!j=i ^j'^ij and x + Y!j=i ^j'^ij ^ -^^x, /?. = 1,...,A;-1. The latter 
relation means that the move Zj^ is applicable to x + ^^.~| e^Zj^. for h = 1, . . . , k. We write 
X ~ y (mod B) if y is accessible from x by B. An accessibility by B is an equivalence 
relation in J^t for any t and each J^t is partitioned into disjoint equivalence classes by B 
(see (fl^ for detail). We call these equivalence classes B-equivalence classes of J-'t. Because 
of symmetry, we also say that x and y are mutually accessible by i3 if x ~ y (mod B). 
Conversely, if x and y are not mutually accessible by B, i.e., x and y are elements from 
two different S-equivalence classes of J?\4x5 we say that a move z = x — y connects these 
two equivalence classes. 

A Markov basis is defined by (j^ as follows. A set of finite moves B = {zi, . . . ,zl} is 
a Markov basis if J-'t itself forms one i3-equivalence class for all t. In other words, if B is 
a Markov basis, every x, y G J?-t are mutually accessible by B for every t. In statistical 
applications, a Markov basis makes it possible to construct a connected Markov chain 
over JF4X for any observed frequency data x. 

Diaconis and Sturmfels (0) showed the existence of a finite Markov basis for any A 
and gave an algorithm to compute one. These results were obtained by showing the fact 
that B = {zi, . . . , zi} is a Markov basis if and only if the set of binomials {u^fc — u^fc , k = 
1, . . . , L} is a generator of the toric ideal associate to A. The algorithm in |3) is based 
on the elimination theory of polynomial ideals and computation of a Grobner basis. 
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2.2 Summary of relevant facts on indispensable moves and min- 
imal Markov bases. 



In ()2l: Il4t Il5fl. we have investigated the question on the minimahty and unique minimahty 
of Markov bases without computing a Grobner basis of J_4. A Markov basis B is minimal 
if no proper subset of i3 is a Markov basis. A minimal Markov basis always exists, because 
from any Markov basis, we can remove redundant elements one by one until none of the 
remaining elements can be removed any further. In defining minimality of Markov basis, 
we have to be careful on signs of moves, because minimal B can contain only one of z or 
— z. Also a minimal Markov basis is unique if all minimal Markov bases coincide except 



for signs of their elements ((Il4l)). 

The structure of the unique minimal Markov basis is given in Here we will 

summarize the main results of the paper without proofs. Two particular sets of moves 
are important. One is the set of moves z with the same value of the sufficient statistic 
t = Az'^ = Az~ , namely 

Bt = {z\ Az^ = Az- = t}, 
and the other is the set of moves with degree less than or equal to n, namely 

Bn = {z I deg(z) < n}. 

Consider the i3|t|_i-equivalence classes of J-'t for each t. We write this equivalence classes 
of J^t as J^t = ^t,i U ■ ■ • U J^t,Kf 

Theorem 2.1 (Theorem 2.1 in (Q)- 

Let B be a minimal Markov basis. Then for 
each t, B n Bt consists of Kt — 1 moves connecting different B\t\-i- equivalence classes of 
J-'t, such that the equivalence classes are connected into a tree by these moves. 

Conversely, choose any Kt — 1 moves zt,i, . . . , zt^ii-^-i connecting different B\t\-i- 
equivalence classes of J-'t such that the equivalence classes are connected into a tree by 
these moves. Then 

[j {Zt,i, . . . ,Zt,i^t_i} 
t:Kt>2 

is a minimal Markov basis. 

From Theorem 12.11 we immediately have a necessarily and sufficient condition for the 
existence of a unique minimal Markov basis. 



Corollary 2.1 (Corollary 2.1 in (Il4l )). A minimal Markov basis is unique if and only 
if for each t, J^t itself forms one B\t\-i- equivalence class or J^t is a two-element fiber. 

This condition is explicitly expressed by indispensable moves. 

Definition 2.1. A move z = z+ — z^ is called indispensable if z^ and z~ form a two- 
element fiber, i.e., the fiber J^ji2+{= J-'j^z-) is written as J-'az+ = {z^,z^}. 
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From the above definition and the structure of a minimal Markov basis, one can show 
that every indispensable move belongs to each Markov basis (Lemma 2.3 in Fur- 
thermore, by the correspondence between moves and binomials, we define an indispensable 
binomial. 

Definition 2.2. A binomial = u^"*" — is indispensable if every system of binomial 
generators of contains or —u^ . 

Clearly, a binomial is indispensable if and only if a move z is indispensable. By def- 
inition, a set of indispensable moves plays an important role to determine the uniqueness 
of a minimal Markov basis: 

Lemma 2.1 (Corollary 2.2 in (0)). The unique minimal Markov basis exists if and 
only if the set of indispensable moves forms a Markov basis. In this case, the set of 
indispensable moves is the unique minimal Markov basis. 

Ohsugi and Hibi further investigated indispensable moves (H;0)- 



Theorem 2.2 (Theorem 2.4 in (j8|)). A binomial is indispensable if and only if 
either or — belongs to the reduced Grobner basis of for any lexicographic term 
order on k[u\ . 

One can find more details in (0). 

3 Definition and some basic properties of indispens- 
able monomials 

In this section we will define indispensable monomials. Then we will show two other 
equivalent conditions for a monomial to be indispensable. We will also prove analogous 
to Theorem 2.4 in (j^, that the set of indispensable monomials is characterized as the 
intersection of monomials in reduced Grobner bases with respect to all lexicographic term 
orders. Hereafter, we say that a Markov basis B contains x if it contains a move z 
containing x by abusing the terminology. 

Firstly we will define an indispensable monomial. 

Definition 3.1. A monomial is indispensable if every system of binomial generators 
of I A contains a binomial f such that is a term of f . 

From this definition, any Markov basis contains all indispensable monomials. There- 
fore the set of indispensable monomials is finite. Note that both terms of an indispensable 
binomial u^^ — are indispensable monomials, but the converse does not hold in general. 

Now we will present an alternative definition. 

Definition 3.2. x is a minimal multi-element if \J-'ak.\ > 2 and |J^^(x-ei)l = 1 foi^ every 
i G supp(x). 
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Theorem 3.1. x is an indispensable monomial if and only if y;. is a minimal multi- 
element. 

Proof. First, we suppose that x is a minimal muhi-element and want to show that it is 
an indispensable monomial. Let n = |x| and consider the fiber JF^^- We claim that {x} 
forms a single ;B„_i-equivalence class. In order to show this, we argue by contradiction. If 
{x} does not form a single S„_i-equivalence class, then there exists a move z with degree 
less than or equal to n — 1, such that 

X + z = (x — z~) + z"*" G J^^x- 

Since |x| = n, \z~\ < n — 1, we have 7^ x — z^ and 

supp(x) n supp(x + z) 7^ 0. 

Therefore we can choose i G supp(x) fl supp(x + z) such that 

^(x — ej) = ^(x + z — ej), X — ej 7^ X + z — ej. 

This shows that \J-'A{x-ei)\ > 2, which contradicts the assumption that x is a minimal 
multi-element. Therefore we have shown that {x} forms a single ;B„_i-equivalence class. 
Since we are assuming that \ J-'ax\ > 2, there exists some other ;B,i_i-class in J-'yu- By The- 
orem |23 each Markov basis has to contain a move connecting a one element equivalence 
class {x} to other equivalence classes of J^^xj which implies that each Markov basis has to 
contain a move z containing x. We now have shown that each minimal multi-element has 
to be contained in each Markov basis, i.e., a minimal multi-element is an indispensable 
monomial. 

Now we will show the converse. It suffices to show that if x is not a minimal multi- 
element, then X is a dispensable monomial. Suppose that x is not a minimal multi-element. 
If X is a 1-element (|J^^x| = 1)) obviously it is dispensable. Hence assume |J?\4x| > 2. In 
the case that is a single S^-i-equivalence class, no move containing x is needed in 
a minimal Markov basis by Theorem 12.11 Therefore we only need to consider the case 
that contains more than one i3n_i-equivalence classes. Because x is not a minimal 
multi-element, there exists some i G supp(x) such that |J^^(x-ei)| ^ 2. Then there exists 
y 7^ X — ej, such that = ^(x — e-j). Noting that |y | = |x — ej| = n — 1, a move of the 
form 

z = y - (x - ei) = (y + e^) - x 
satisfies < deg(z) < n — 1. Then 

y + ei = X + z 

and X and y + e-j belong to the same ;B„_i-equivalence class of J\4x- Since x 7^ y + e^. 
Theorem 12. II states that we can construct a minimal Markov basis containing y + e^, but 
not containing x. Therefore x is a dispensable monomial. □ 

We will give yet another definition. 
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Definition 3.3. ^ is a minimal i-lacking 1-element if \J-'ax\ = 1; |-^^(x+ei)| > 2 and 
|-^^(x+e,-e,)| = 1 for each] G supp(x). 

We then have the following result. 

Theorem 3.2. The following three conditions are equivalent 1) ^ is an indispensable 
monomial, 2) for each i G supp(x), x — ej is a minimal i-lacking 1-element, 3) for some 
i e supp(x), X — ej zs a minimal i-lacking 1-element. 

By the previous theorem we can replace the condition 1) by the condition that x is a 
minimal multi-element. 

Proof. 1) =^ 2). Suppose that x is a minimal multi-element. Then for any i G supp(x), 
X — ej is a 1-element and |JF4(-(x_ei)+ei)| = 1-^^x1 > 2. If x — e^ is not a minimal i-lacking 
1-element, then for some j G supp(x — ej), |J^^(x-ej)| > 2. However j G supp(x — e^) C 
supp(x) and \ J-'A{x-ej)\ > 2 contradicts the assumption that x is a minimal multi-element. 
It is obvious that 2) =^ 3). 

Finally we will prove 3) =^ 1). Suppose that for some i G supp(x), x — e^ is a minimal 
i-lacking 1-element. Note that \J-'ax.\ = \^Ai{y.-ei)+ei)\ > 2. Now consider j G supp(x). 
If j G supp(x - ej) then |JU(x-e,)| = \J^A{(^-e,)+e,-ej)\ = 1- On the other hand if 
j ^ supp(x — ej), then j = i because j G supp(x). In this case \ J-'A{y.-ei) \ = 1- This shows 
that X is a minimal multi-element. □ 

Theorem 13.21 suggests the following: Find any 1-element x. It is often the case that 
each ej, z = 1, . . . ,p, is a 1-element. Randomly choose 1 < i < p and check whether x-|-ej 
remains to be a 1-element. Once |J^x+eJ > 2, then subtract e^'s, j ^ i, one by one from 
X such that it becomes a minimal z- lacking 1-element. We can apply this procedure to 
finding indispensable monomials of some actual statistical problem. 

For the rest of this section we will illustrate this procedure with an example of a 
2x2x2 contingency table. Consider the following problem where p = 8, d = 4 and A is 
given as 

/ 1 1 1 1 1 1 1 1 \ 
. _ 1 1 1 1 
110 110 
\ 1 1 1 1 / 

In statistics this is known as the complete independence model of 2 x 2 x 2 contingency 
tables. To see the direct product structure of I explicitly, we write indeterminates as 

U = (Mill, Mll2, Ml21, Ml22, ^211, ^212, ^221, ^222)- 

To find indispensable monomials for this problem, we start with the monomial = um 
and consider x + ej , z G X. Then we see that 

• Uiii, niiiMii2, ^111^121, M111U211 are 1-elements, 

• ^111^122, W111U212, W111M221 are 2-elements and 
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• ^111^222 is a 4-element. 



From this, we found four indispensable monomials, 'Uiii'Ui22, M111M212) Wiii-U22i and U111U2221 
since each of M122, '"212, ^221, M222 is a 1-element. 

Starting from the other monomials, similarly, we can find the following list of indis- 
pensable monomials, 

• U111U122, Mlll^i212, M1IIM22I, ^ill2Ml21, Mll2^i211, M112M222, ^il2lM222, Ml21^i211, 

'Ui22'U22i, 'Ui22'U2i2, M211M222, ^212^221, each of which is a 2-element monomial, and 

• ^^111^222, Mii2^i22i, M121M212, ^^1221*211, each of which is a 4-element monomial. 

The next step is to consider the newly produced 1-element monomials, 
'^111, M111M112, 'itiii'Ui2i, 'Uiii'it2ii and so on. For each of these monomials, consider adding 
ej,i G X one by one, checking whether they are multi-element or not. For example, we 
see that the monomials such as 

«iii'«iii^ii2,^iii«n2' • • • 

are again 1-element monomials (and we have to consider these 1-element monomials in 
the next step). On the other hand, monomials such as 

uluUi22, UiuUu2Ui22, ^111^222, ^111^^ 112^221 , • • • 

are multi-element monomials. However, it is seen that they are not minimal multi-element, 
since 

^111^^122, M112M122, ^^111^222, ^^112^^221, ■ ■ ■ 

are not 1-element monomials. 

To find all indispensable monomials for this problem, we have to repeat the above 
procedure for monomials of degree 4, 5, . . .. Note that this procedure never stops since 
there are infinite 1-element monomials, such as 

"ill "1125 • • • 



for arbitrary n,m. This is the same difficulty mentioned in Section 2.2 in (|lj). Since 
indispensable monomials belong to any Markov basis, in particular to the Graver basis. 
Theorem 4.7 in (fl3) gives an upper bound for the degree of indispensable monomials and 
we can stop at this bound. 

Finally we will prove the following theorem, which is analogous to Theorem 2.4 in 
but much easier to prove, since it focuses on a single monomial (rather than a binomial). 
We need to reproduce only a part of the proof for Theorem 2.4 in (j^. 

Theorem 3.3. A monomial x is indispensable if for every lexicographic order <iex the 
reduced Grobner basis with respect to <iex contains x. 
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Proof. It suffices to show that if a monomial x is dispensable, then there exists a lex- 
icographic term order <iex such that the reduced Grobner basis i3<j^^ does not contain 
x. Note that the positive part and negative part of a move belong to the same fiber. 
Therefore if x is a 1-element, then no Markov basis contains x. In particular no Grobner 
basis contains x. Therefore we only need to consider x such that |jF4x| > 2. 

Since we are assuming that x is dispensable, there exists a Markov basis B, which 
does not contain x. Then there exists a move z = z+ — G i3, with z's sign changed if 
necessary, such that z is applicable to x, i.e., x > z . Since B does not contain x, z" 7^ x 
and hence z~ is strictly smaller than x. Now choose <iex such that the initial term of z 
is z~. Then the reduced Grobner basis with respect to <iex does not contain x. □ 



4 Further properties of indispensable monomials 

In the previous section we gave some basic characterizations of indispensable monomials. 
In this section we will show further properties of indispensable moves in terms of minimal 
Markov bases in (14) and a norm-reducing Markov basis in ([isl). 

Firstly, we will state the following lemma, which is already implicitly used in the proof 
of Theorem 13.11 

Lemma 4.1. A monomial is indispensable if and only if J^ax contains more than one 
B\^\^i- equivalence class and the one-element set {x} forms a B\^\^i - equivalence class. 

Proof. Suppose that JF^^ contains more than one i3|x|-i-equivalence class and the one- 



element set {x} forms a i3|x|-i-equivalence class. Then by Theorem 2.1 in (|lj), every 
Markov basis has to connect x with some other i3|x|-i-equivalence class of J-'ax- Therefore 
X has to appear as a positive part or a negative part of some move z of the Markov basis. 

Conversely, we show that if J^^x contains just one i3|x|-i-equivalence class or the 
equivalence class containing x contains some other vector y, then x is dispensable. In the 
former case, JF4X is already connected by moves of degree less than or equal to |x| — 1 and 
no minimal Markov basis contains a move having x as the positive or the negative part. 
On the other hand if y 7^ x belongs to the same S[x|-i-equivalence class, then by Theorem 



2.1 in ()l4). there exists a minimal Markov basis involving y and not x. Therefore x is 



dispensable. □ 



From ()l4 ) it follows that the moves of all minimal Markov bases bel ong to a common 



=et of flbe~Also, we defined the minimum fiber Markov ba.es Bmp in P) as 

Bmf = {z = z+ - z" I z+ 7^ z" (mod i3|z|-i)}- 

Based on Lemma 14.11 now we will prove four propositions concerning the fibers in i^MP- 
In the following four propositions, an equivalence class of a fiber J-'t means a B\t\-i- 
equivalence class of J^f 

Proposition 4.1. The following three conditions are equivalent: 1) all equivalence classes 
of all fibers of Bmf o,re singletons, 2) there exists a minimal Markov basis, such that all 
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monomials contained in the basis are indispensable. 3) for all minimal Markov bases, all 
monomials contained in the basis are indispensable. 

Proof. Obviously 3) ^ 2). 2) ^ 1) follows from Lemma f4. II because a minimal basis has 
to connect all equivalence classes of each fiber of Bmf into a tree. To show that 1) ^3), 
we again use the fact that a minimal basis has to connect all equivalence classes of each 
fiber of i^MF into a tree. If all equivalence classes of a fiber are singletons, then both 
terms of a move connecting two equivalence classes are indispensable. This completes the 
proof. □ 

Proposition 4.2. There exists a minimal Markov basis such that each move of the basis 
contains an indispensable monomial if and only if each fiber of Bmf contains a singleton 
equivalence class. 

Proof. Let ;B be a minimal Markov basis such that each move of B contains an indispens- 
able monomial. This monomial forms a singleton equivalence class. Therefore each fiber 
of Bmf contains a singleton equivalence class. Conversely if each fiber of Bmf contains a 
singleton equivalence class, we can construct a tree which connects each equivalence class 
of the fiber to the singleton equivalence class. Then the resulting minimal Markov basis 
has the property that each move of the basis contains an indispensable monomial. □ 

Proposition 4.3. Every move of any minimal Markov basis contains an indispensable 
monomial if and only if all but one equivalence classes of each fiber of Bmf o'^e singletons. 

Proof. If all but one equivalence classes of each fiber of Bmf are singletons, then in con- 
necting these equivalence classes into a tree, each move has to contain an indispensable 
monomial. On the other hand if there exist two non-singleton equivalence classes in a 
fiber, then we can construct a minimal Markov basis containing a move connecting these 
two equivalence classes. This move does not contain an indispensable monomial. □ 

Next, we consider indispensable monomials in terms of norm-reduction introduced in 
(jlih . We will give a definition of a norm-reducing Markov basis here (see fllih for detail). 

Definition 4.1. A set of moves B is 1-norm reducing if for all t and for all x, y G J?-t 

with X 7^ y, there exist some z & B and e G { — 1, +1} satisfying either 

|(x + ez) -y| < |x-y| 

or 

|x- (y + ez)| < |x-y|. 

It is easy to show that, if B is 1-norm reducing, then it is a Markov basis (see Propo- 
sition 1 in (jlil )). Therefore we call B a 1-norm reducing Markov basis if it is 1-norm 
reducing. An example of 1-norm reducing Markov basis is the Graver basis (see Proposi- 
tion 2 in (fish). Now we will give a characterization of indispensable monomials in terms 
of the norm reduction. 
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Proposition 4.4. A move, whose both monomials are indispensable, belongs to each 1- 
norm reducing Markov basis. 

Proof. Let z = z+ — be a move such that both z+ and z~ are indispensable, i.e., 
{z+} and {z~} are singleton equivalence classes of a fiber. If a Markov basis E does not 
contain z, then it can not decrease the distance between z"*" and z^. Therefore B is not 
1-norm-reducing. □ 

Finally, we give a further definition, which is similar to a minimal multi-element. 

Definition 4.2. Tt is a minimal multi-element fiber if { \Tt\ > 2 and |t| = 1 } or { 
l-^tl > 2, |jFtJ = iJFt^l = 1 for any t = ti + 12 satisfying |J-tJ, \Tt^\ > 1 }■ 

The meaning of this definition is as follows. Suppose \ J-'t^ \ > 2 and xi, X2 G J^ti- Then 
for any t2 = ^X3, |jFtj_,_t2| > 2 follows since Xi + Xa, X2-I-X3 G J^-t^+tj. Note that the former 
case, \T\ > 2, |t| = 1, corresponds to the case Tt 3 e,, e^, . . ., for some z, j, . . .. One of the 
situations that this special case appears is ^ = [!,...,!]. Hereafter, we only consider the 
problem that |t| > 2 holds for every minimal multi-element fiber J^m- In other words, 
we assume that every is a one-element. In this case, minimal multi-element fiber is 
characterized as follows. 

Proposition 4.5. J^t is a minimal multi-element fiber if and only if all the elements in 
Tt are indispensable monomials. 

Proof. Suppose all the elements in J-'t are indispensable monomials and t = ti + 12 where 
l-^til > 2, iJFtJ > 1. Write Xi,X2 G J-'t^ and X3 G J-'t2- In this case, both Xi + X3 and 
X2 + X3 are in J-'t and therefore indispensable monomials by the assumption. However, 
for any i G supp(x3), xi -|- X3 — ej and X2 -|- X3 — are again in the same fiber, which 
contradicts the assumption that xi + X3 and X2 + X3 are minimal multi-elements. 

Conversely, suppose J-'t is a minimal multi-element fiber and x G is dispensable. 
In this case, since x is not a minimal multi-element, there exists some satisfying 
\J^A{x-ei)\ ^ 2. Therefore we have ^x = ^(x — Cj) -|- Aei, i.e., t = ti + t2 where 
ti = ^(x — Bj) and t2 = Aei, which contradicts the assumption that J^t is a minimal 
multi-element fiber. 

□ 

5 Examples 

In this section, we will give some indispensable monomials and dispensable monomials in 
minimal Markov bases for some statistical models. As is stated in fjl4). there are some 
models where a minimal Markov basis is uniquely determined, and some models where it is 
not uniquely determined. Furthermore, by considering the indispensability of monomials 
contained in minimal Markov bases, we can classify Markov bases by the indispensability 
of monomials as follows. 
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• Case 1. A minimal Markov basis is uniquely determined, i.e., the set of indispensable 
moves forms a Markov basis. 

• Case 2. A minimal Markov basis is not uniquely determined, but all the monomials in 
minimal Markov bases are the same and indispensable. In this case, all equivalence 
classes of each fiber of Bmf are singletons. 

• Case 3. A minimal Markov basis is not uniquely determined, and they contain some 
moves where their positive or negative parts are dispensable monomials. In this 
case, some equivalence classes of some fiber of i^MF are not singletons. 

We will show examples for Case 2 and Case 3 in this section. As for Case 1, the set 
of the positive and negative parts of indispensable binomials is the set of indispensable 
monomials. One of the most simple examples for Case 1 is an independence model of two- 
way contingency tables. A quite difficult example is a no three-factor interaction model 
of three-way contingency tables, i.e., the case that is the two-dimensional marginal 
totals of three-way contingency tables x. For this example, minimal Markov bases for 
some small sizes of x is shown to be unique (see (0) for example of 3 x 3 x case). 
Indispensable monomials for Case 1 clearly coincide the positive and negative parts of 
indispensable binomials. 



5.1 Examples of Case 2 

One-way contingency tables with fixed totals. First we consider the simplest 
example given by ^ = lp,p > 2, where Ip = (1,...,!)' is the p dimensional vector 
consisting I's. As is shown in (|1J), minimal Markov bases for this problem contain 
dispensable moves only, which connect p elements, 

{Mi,M2, ■■■,Up} 

into a tree. It is also obvious that these p monomials are all indispensable. 



Complete independence models of three-way contingency tables. 

a generalization of the problem considered at the end of Section El 

Let X be a frequency vector for I x J x K contingency tables and let 

X = {ijk I 1 <i < 1,1 < j < J,l < k < K}. 



We will show 



A is given as 



A 



I'j ® E_ 



K 



K 



K 



where En is the nxn identity matrix. The minimum fiber Markov basis for this problem 
is given in ([ish as 
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BlBP = {Uijiki 

U \UiijikV'i2j2k '^i\j2k'^i2jiki ^1 7^ ^2; Jl 7^ J2}) 
^ {_^i\jik\'^i2j2k2 ^njifc2^i2j2fci ) '^iij\k\'^i2j2k2 '^iij2k\'^'i2jik2i 

^iljlfci '^12^2^2 '^iij2k2^i2jlkij 'I^iijik2'^i2j2ki '^i\j2k\'^i2jik2i 
'^hjlk2'^i2j2kl ^ili2fc2^j2ilfcl 5 ^ni2fcl^*2jlfc2 ^41^2^2 ^*2ilfel 5 

^1 ^ «2, jl 7^ J2,fci 7^ A;2}- 

Here, i3iDP is the set of indispensable moves. i3* is the set of all degree 2 moves which 
connect all elements of the four-elements fiber 

ili2jij2kik2 {-^ {"^ufc} I "^ii-- •^i2-- '^■ji' '^■j2- •^■■ki ■^■■k2 

{'^iljlkl'^i2j2k2^ '^iljlk2'^'i2j2kl1 '^11^2^1 ^*2ilfe2 J ^nj2fc2'^*2ilfel }• 

The minimal Markov basis in this case consists of Bijy^ and three moves for each ii 7^ 
^2,^1 7^ ^2, and ki 7^ which connect four elements of ^ni2jij2kik2 iiito a tree. In this case, 
the four elements of ^iii2jij2kik2 are different i3i -equivalence classes, which are obviously 
singletons. Therefore the set of indispensable monomials for this problem is 

{'^iijikiUi2j2k2J '^iijifc2^j2i2fen '^i\j2ki'^i2jik2i '^iij2k2^i2jlk-ii ^1 7^ ^2; jl 7^ j2) ^1 7^ ^2} 

in addition to the positive and negative components of i3iDp. Figure ^ illustrates the fiber 

•^iii2jij2kik2 ■ 




Figure 1: Illustration of the 4-element fiber ^iii2jij2kik2 of the complete independence 
model of three-way contingency tables. The four circles are four monomials, where xi = 
^*iii'si^*2j2's2 5 -^^2 ^iijik2^i2j2ki}^3 ^iij2ki^i2jik2} and X4 ^iij2k2^i2jiki- Each mouomial 
forms i3i-equivalence class of the fiber by itself represented by the dotted square. We 
will use this convention in all forthcoming figures. The thick lines mean a choice of 
three dispensable moves, {xi — X2,xi — X4,X3 — X4}, which is an example of choices for 
constructing a minimal Markov basis. 
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Hardy- Weinberg model. Another example considered in ()l4^ is the Hardy- Weinberg 
model for / alleles, i.e., 



X — (a^ii, a;i2, . . . , xu, X22, . . . , X2/, 0:33, . . . , xu)' 



and 



A = {Ai Ai-i ■ ■ ■ Ai), Ak = {Okx{i-k) -Bfc) : 
where is the following k x k square matrix 



Bk 



2 1 1 ■■■ 1 
1 ■•■ 

1 











As is stated in (|l4). a minimal Markov basis for this case is not unique, and the minimum 
fiber Markov basis is: 

Bmf = i3iDP U B*, 

^IDP {'^iiii'^i2i3 '^'i\i2'^iiizi '^i\il^i2i2 ^1-112^1 

where ii,i2, is, h are all distinct, and Uij = Uji for i > j. Here, B* is the set of all degree 
2 moves which connect all of the elements of the three-element fiber 



III 



{^iii2^i3i4,, ^iiiz^i2iii ^1114^1213}' 



Again, these three elements of Ti^i^i-j^ii form singleton ;Bi-equivalence classes of it by 
themselves, and are indispensable monomials. Figure |21 illustrates the fiber J^ij^i2i3i4- 



5.2 Examples of Case 3 

Some examples for this case are found in the hierarchical models of 2x2x2x2 contingency 
tables considered in (jl|). First we will show one of them as an example of Case 3. By 
modifying the example, we will show another example of the situation considered in 
Proposition l4.3l i.e., the situation that some dispensable moves contain both indispensable 
and dispensable monomials as their positive and negative parts. 



12/13/23/34 model of 2x2x2x2 contingency tables. Let x be a frequency vector 
for 2 X 2 X 2 X 2 contingency tables (p = 16). We write indeterminates with respect to a 
lexicographic order as 

U = {Wllll, M1112, M1121, ^1122, ^1211, • • • , ^^2222}- 
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Figure 2: Illustration of the 3-element fiber Tiyi^i^i^ of the Hardy- Weinberg models. The 
three circles are three monomials, Xi = -Uj^jj-Uigj^, X2 = Ui^i^Ui^i^^yi.j, = Ui^i^Ui^i^. The thick 
lines are the two dispensable moves, {xi — X2, X2 — X3}. 



Consider the model of d = 9 given as 



A 



1111111111111111 
1111111100000000 

1111000011110000 
1100110011001100 





1 
1 



1 

1 1 





10 10 



1 



10 10 1 
000000000000 
110000000000 
1100000011000000 
1000100010001000 



For this case, there are 12 indispensable moves of degree 2 and 4 indispensable moves of 
degree 4, but the set of indispensable moves does not form a Markov basis. In addition 
to the indispensable moves, we have to consider moves connecting 4-element fiber 

— {'iillll'iil22l'ii212l1i2212, 1illl2'iil22l'ii212l1i2211 , 1ill2l1il21l1i2112'ii2221 , 1*11211*12121*21111*2221} 

and 8-element fiber 

= {ltllll1tl22l1t2122lt2212 5 1*11121* 12221*21211*2211 , 1*1111^12221*21211*2212 ; 1*11121* 122l1*2122l*2211 ; 

1*112l1*121l1*2112l*2222, 1*1122l*1212l*211l1*2221 , 1*112l1*1212l*211l1*2222 , 1*1122l*121l1*2112l*222l}- 

For JF^ , we have vBs-equivalence classes of it as 

— {l*llll1*122l1*212l1*2212, 1*1112l*122l1*212l1*2211 }U{Wll2l1*121l1*2112l*2221 , 1*112l1*1212l*211l1*222l}- 

Therefore these 4 elements are dispensable monomials. In fact, we can find an element of 
minimal Markov basis not containing ^1111^12211*21211*2212, for example, as 



1*1112l*122l1*212l1*2211 ~ 1*112l1*121l1*2112l*2221 • 
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Similarly for T'^ , yBs-equivalence classes of it are given as 

= {WllllMl22l'W2122'W2212; 'Wlll2'Wl222'W2121 "^2211 ; 'Wllll'Wl222'W212l'W2212 j Wiii2ltl22l1t2122't*2211 } ; 

U{nii21'Ul21l'W2112'U22225 1*1122^^1212^2111^^2221 , Mii2lMl212W2111^i2222 , 1^1122'Wl21l1l2112W2221 } • 

Figure 01 and Figure lU illustrate the fiber JF^ and JF^, respectively. 
















f X4) 



Figure 3: Illustration of the 4-element fiber of 12/13/23/34 model of 2 x 2 x 
2x2 contingency tables. The four monomials are Xi = MiiiiMi22iM2121'W2212, X2 = 

Mlll2Ml22lM212lM2211,X3 = Mll2lMl21lM2112M2221 , and X4 = Mll2lMl212M211lM2221 • {xi,X2} and 

{x3,X4} form iSa-equivalence classes of the fiber. The thick line is a dispensable move, 

Xi - X4. 



12/13/23/34 model of 2x2x2x2 contingency tables with a structural zero cell. 

We modify the previous example by introducing a structural zero cell, xm = 0. This 
situation corresponds to removing the indeterminate Uim and the first column of A as 



u 



{Miii2, "^11217 ^1122) W1211, . . . , U2222\i 



A 



111111111111111 
111111100000000 
000011110000 
110011001100 
101010101010 
111000000000000 
100110000000000 
100000011000000 



1 1 1 

1 
1 



000100010001000 
In this case, the fiber in the previous example is modified to a 3-element fiber. 



T * = {Mlll2Ml22l1i212l1i2211,1ill2lMl21lM2112?i2221,1^112l1il212M211lM222l}- 



17 








® 


® 




® 


® 


® 





V. > 



Figure 4: Illustration of the 8-element fiber J^^ 12/13/23/34 model of 2 x 2 x 
2x2 contingency tables. The eight monomials are xi = ^^1111^^1221^2122^22125 

X2 = Miii2Mi222'W2121^2211; X3 = Mllllttl222^i212lW2212; X4 = Ulll2Wi22lW2122^2211 ; X5 = 
^1121^^1211^2112^2222; = Mll22^1212'W2111^2221 ; X7 = Mll2lMi212^2111^2222 5 ^nd Xg = 

Mii22^ii2ii«2ii2«222i- {xi, X2, X3, X4} and {X5, Xg, X7, Xg} form Ss-equivalence classes of 
the fiber. The thick line is a dispensable move, X2 — xg. 



Since JF has still different Ss-equivalence classes, we have to consider moves connecting 
the elements of J-'^* to construct a minimal Markov basis. In this case, i^a-equivalence 
classes of JF^* are given as 

= {^illl2^^1221^i2121^i221l} U {^1121^1211^2112^2221 , ttll21^il212^^211lM2221 } • 

Therefore a minimal Markov basis for this problem has to contain either 

'?illl2^il221'lf212l'lt2211 — 1fll2lMl21lM2112^i2221 

or 

'Ulll2^1221^212l'it2211 " ^1121^1212'it211lM2221 ■ 

The above two moves are dispensable, and the negative parts of both moves are also 
dispensable monomials, whereas the positive part, 'Uiii2'Ui22iW2i2i'W22ii, is an indispensable 
monomial. Figure El illustrates the fiber JF^*. 

6 Some discussions 

In this paper, the concept of indispensable monomials is introduced, by extending the 
notion of indispensable binomials. Both in the framework of Markov bases and toric 
ideals, the indispensable monomial plays an important role since it has to be included in 
all Markov bases or generators of toric ideals. It is true that enumerating indispensable 
monomials is as difficult as enumerating indispensable binomials. 

Note that, by the notion of indispensable monomials, we can characterize a dispensable 
binomial as (i) a difference of two dispensable monomials, (ii) a difference of dispensable 
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Figure 5: Illustration of the 3-element fiber JF^* of 12/13/23/34 model of 2x2x2x2 
contingency tables with structural zero cell. This fiber is constructed by removing element 
xi from . In this case, {X2} and {x3,X4} form i^s-equivalence classes of the fiber. 
The thick line means a choice of a dispensable move, X2 — X4. Another possibility of 
constructing a minimal Markov basis is to choose a dispensable move, X2 — X3. The 
monomial X2 is included in any minimal Markov basis, and is an indispensable monomial. 



and indispensable monomials, or (iii) a difference of two indispensable monomials. The 
situations where each case arises are shown in Proposition 14.31 and in Proposition 14.11 We 
have found some examples for the case (ii) by introducing some structural zero cells for 
the case (i). 

The enumeration of indispensable monomials seems very important problem, since it 
can lead directly to the enumeration of indispensable binomials. In addition, it also gives 
the fibers of the special structure that it contains at least one singleton equivalence class. 
Moreover, by finding dispensable binomials which are differences of two indispensable 
monomials, we can find all fibers that only contain singleton equivalence classes. 
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